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Abstract 

We study the behaviour of the magnetization in a half-metalhc ferromagnet /nonmagnetic insula- 
tor/ferromagnetic metal/paramagnetic metal (FM1/NI/FM2/PM) tunnel junction. It is calculated 
self-consistently within the nonequilibrium Keldysh formalism. The magnetic regions are treated as 
band ferromagnets and are described by the single-band Hubbard model. We developed a nonequi- 
librium spectral density approach to solve the Hubbard model approximately in the switching 
magnet. By applying a voltage to the junction it is possible to switch between antiparallel (AP) 
and parallel (P) alignment of the magnetizations of the two ferromagnets. The transition from AP 
to P occurs for positive voltages while the inverse transition from P to AP can be induced by neg- 
ative voltages only. This behaviour is in agreement with the Slonczewski model of current-induced 
switching and appears self-consistently within the model, i.e. without using half-classical methods 
like the Landau-Lifshitz-Gilbert equation. 
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There has been considerable interest in the phenomenon of current-induced switching of 
magnetization since it was first proposed over 10 years ago^i^. The basic idea behind this 
effect is as follows: the spin direction of electrons moving in a ferromagnet will be mostly 
aligned parallel to the magnetization axis. When these spin-polarized electrons are trans- 
ported to a second ferromagnet, e.g. by applying a voltage, then the spin angular momentum 
of the itinerant electrons will exert a torque on the local magnetic moment. This torque 
is known as the spin-transfer torque. It will have an influence on the direction of magne- 
tization. If the parameters of the materials are chosen in the right way and if the current 
through the junction is high enough it is even able to switch the magnetization of one ferro- 
magnet from parallel to antiparallel or vice versa relative to the other one. This effect was 
seen both in all-metallic junction^i^s^ such as Co/Cu/Co and in magnetic tunnel junctions 
(MTJs) consisting of two ferromagnets divided by a thin nonmagnetic insulator-i^i^. In this 
paper we focus on a special case of the latter, where the ferromagnetic lead is half-metallic, 
i.e. there are only electrons of one spin direction present at the Fermi energy. 
Some of the possible technological applications of spin-transfer torques in MTJs have been 
discussed by Diao et al.^. Most of the theoretical work in this area of research has been 
focused on the Landau-Lifshitz-Gilbert (LLG) equationi^'ii'i^'i^'i^, which is a macroscopic, 
half-classical equation. The torques entering this equation were usually calculated in a mi- 
croscopic picture while treating the interactions on a mean field level. In this paper we 
propose a model which takes interactions beyond mean field into account. We make no use 
of the LLG equation or other macroscopic approaches and thus we stay on the quantum 
mechanical level throughout this paper. 

We will start the presentation of the theory by introducing a model Hamiltonian which 
describes the magnetic tunnel junction shown schematically in Fig. [H There are two fer- 
romagnetic metals {L and R) divided by a non-magnetic insulator (/) and additionally a 
paramagnetic metal (P) which is necessary to have a well-defined chemical potential on the 
right side of the second ferromagnet. Each region consists of a single s-like band. The two 
outer leads L and P are treated as semi-infinite. 
The total Hamiltonian consists of several parts: 

H = Hl + Hu + Hj + Hri + Hn + Hrp + Hp (1) 

Hl{r) describes the left (right) ferromagnet, Hj the insulator and Hp the paramagnet. Both 
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FIG. 1: Schematic picture of the magnetic tunnel junction with apphed voltage V. The conduction 
bands are shown as rectangles. Occupied states in the metals are hatched and the directions of 
magnetization in the ferromagnets are symbolized by thick arrows. 

insulator and paramagnet are assumed to be non-interacting, so their Hamiltonians consist 
of the kinetic energy only: 



c^kxo-(c^kxo-) annihilation (creation) operator of an electron with wave vector kx and 

spin a. ekj^ is the dispersion of the lattice which throughout this paper is chosen as a tight- 
binding bcc lattice. The applied voltage V will shift the center of gravity of the paramagnet 
by Vp = \/ and half of that amount for the insulator, V/ = V jl. Positive voltage, V > 0, will 
shift the bands to lower energies, while negative applied voltages result in a shift towards 
higher energies. 

The Hamiltonians of the left (L) and right (R) ferromagnet are formally almost identical. 
Besides the kinetic energy they also include on-site Coulomb interaction. They are given in 
a mixed Bloch-Wannier representation: 



where M stands for either L or R. The Hubbard-C/ determines the interaction strength. 
^iMO- — ^1^M<y^'i-MC is the occupation number operator. The voltage V shifts only the band 
center of the right ferromagnet, Vr — V-, while the left ferromagnet is not directly influenced 




(2) 




(3) 
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by V, i.e. Vl = 0. 

The remaining three terms of the Hamiltonian are responsible for the couphng between the 
different regions. These couphngs act as a hybridization^^ between the bands and therefore 
the Hamiltonians are (M = L, R; X = I , P): 

Hmx = Yl {^^Mkx4Mad\.xa + h.c.) (4) 

kMkxo- 

They are characterized by the couphng constants Ck^^kx which determine the strength of 
the hybridization between the different bands. In general the couplings are wave vector 
dependent but for the sake of simplicity we neglect this dependence, ei^^kx = ^mx = ^xm- 
Furthermore we assume the coupling between the ferromagnets and the insulator to be equal, 
so that eu = ^ri = ^mi- Altogether there remain two couplings 64/7 and eRp which cannot 
be calculated within this model, so they will be treated as parameters. 
The main topic of this work will be the calculation of the non-equilibrium magnetization 
m of the right ferromagnet within the Keldysh formalism!^. It can be calculated with the 
help of the Fourier transform of the so-called lesser Green's function defined as G^^^{t, t') = 



m = — n\ = 



^ / dEY,{G^AE)-G^,,{E)) (5) 

-00 

rifj = (ficr) is the occupation number of particles with spin a in the right ferromagnet. 
In order to derive the lesser Green's function one first has to calculate the retarded one, 
^k_BO-(-^) — ((ckjjo-; c-kRcr))^- By using the equation of motion method one finds 

Two different selfenergies appear in the Green's function. First there is the interaction 
selfenergy which can only be calculated approximately for the Hubbard model. We propose 
a non-equilibrium spectral density approach (NSDA). The basic idea behind this approach is 
to choose the selfenergy in such a way that the first four spectral moments are reproduced by 
the theory. Some details of its derivation are given in Appendix |Xl The mean field (Stoner) 
solution of the Hubbard model on the other hand satisfies only the first two moments. One 
finds for the selfenergy 
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This expression is coincidentally formally identical to the equilibrium spectral density 
approach^. The difference is in the spin-dependent band correction B^cr which is given 
by 



^ -oo 

+ - l) - To,r) [E - eu, - T^l^-AE)) ]g<^_AE) + ^A<^_(i?)j8; 

It has to be calculated self-consistently since also appears on the right hand side as 
part of the lesser Green's function. Tq ^ is the center of gravity of the right ferromagnet. 
The second selfenergy is the transport selfenergy which is due to electrons hopping between 
the different materials. Its retarded (lesser) component is given by 



kj kp 



G\^^f{E) is the Green's function of the insulator when it is only coupled to the left ferro- 
magnet, i.e. 

Gl'^fiE) = — ' . . (10) 

Since we neglected the wave vector dependence of the couplings, the transport selfenergy 
is only formally dependent on the wave vector. glf_j^„{E) and gl^pAE) are the equilibrium 
Green's functions of the left ferromagnet and the paramagnet, respectively. They can be 
easily calculated by the equation of motion method. One finds for M = L,P: 

The paramagnet does not include interactions so that SJ^^^ = —iO^. Since we are mainly 
interested in the properties of the right ferromagnet, we assume that the left one is half- 
metallic so that its minority states play no role for small voltages. This is done by using the 
mean field selfenergy SJ^^^ = UlUl-ct with sufficiently large Ul- Thus the retarded Green's 
function is known. 

The lesser Green's function follows immediately from the Keldysh equation: 

G^.AE) = Gl^AE)A<^AE)Gl^AE) (12) 
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where the advanced Green's function is simply the complex conjugated of the retarded 
one, G^^^{E) = {Gl^^^{E))* . Furthermore we need the lesser component of the transport 
selfenergy which was already defined in Eq. The lesser part of the insulator Green's 

function can again be calculated with the help of the Keldysh equation: 

GiTiE) = ^ E Gi'^l\E)el,g<^^{E)Gi'^f{E) (13) 

Since the Green's functions in the left ferromagnet and the paramagnet are equilibrium 
quantities, their lesser parts read 

9l,,,AE) = -2zh^p){E)lmgl^^^^^{E) (14) 

where fL{p){E) is the Fermi function in lead L{P) with chemical potential ^l{p)- They 
are related by Hl — l^'P = V. Thus we have a closed set of equations for calculating the 
magnetization of the ferromagnet. 

In Fig. [2] a typical numerical solution for the voltage-dependent magnetization is shown. 
We will first discuss the black curve, which was calculated with a hybridization strength of 
^Mi = 0.5 eV. For the calculation we started with parallel alignment of the two magnetiza- 
tions (point A in the figure). Then a negative voltage is applied, i.e. the right ferromagnet 
is shifted to higher energies compared to the left one. At a critical voltage the parallel 
alignment becomes unstable and the magnetization reverses its sign {B to G). Thus the 
magnetizations are now antiparallel. When the voltage is further decreased the magneti- 
zation stays more or less constant until point D is reached. Then the process is reversed 
and the voltage is reduced to zero again. The magnetization follows the same line as before 
until the switching point G is reached. There it does not switch back to parallel alignment 
but rather stays at about the same level. When E is reached the direction of the voltage 
is reversed, i.e. the right ferromagnet will now be shifted to lower energies. For small volt- 
ages there is only a slight increase until a critical voltage is reached (F). This voltage has 
approximately the same value as the first one at point B, but of course with an opposite 
sign. There the antiparallel alignment is no longer stable and the systems returns to its 
initial parallel state which is not influenced by higher voltages {G to H). Then the voltage 
is turned off and the system will be at its starting point A again, so the hysteresis loop is 
complete. 

As another test we start again at point A, but this time we turn on a positive voltage. Then 
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FIG. 2: Numerical results for the magnetization as a function of applied voltage for two different 
values of the hybridization strength emi between the metals and the insulator. The arrows indicate 
in which direction the voltage was changed. Parameters: Band occupation n = 0.7; band widths: 
Wl = 3 eV, Wi = l eV, Wr = 2 eV and Wp = 5 eV; interaction strengths: Ur = 4 eV, C/l = 4 eV; 
band center of the insulator Tqj = 5 eV; hybridization strength: €rp = 0.05 eV; temperature: 
T = K 

no switching occurs and the system will move reversibly to point H. A similar reversible 
behaviour is seen when the alignment is antiparallel {E) and the voltage is decreased. This 
is shown by the arrows in the figure. So, one has to conclude that switching of magnetization 
from parallel to antiparallel alignment is only possible for negative voltages and the reverse 
process will only appear for positive voltage. The behaviour just described is one of the 
hallmarks of current-induced switching of magnetization and thus our proposed model is 
indeed able to simulate this effect without leaving the microscopic picture! 

Now we want to give a short explanation on how exactly our model is able to provide 
these results. The key to the understanding lies in the effect the hybridization parts of the 
Hamiltonian have on the quasiparticle density of states (QDOS) of the switching magnet 
and in the polarization of the current. A hybridization between two bands generally will 
lead to a repulsion between them, i.e. the energetic distance between their respective centers 
of gravity will increase the stronger the effect of the hybridization is. The magnitude of this 
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FIG. 3: (Color online) Quasiparticle density of states for the parameter set of Fig. [2] for F = 
and eMi = 0.5 eV. Upper picture for parallel alignment of the magnetizations (point A), lower 
picture for antiparallel alignment (point E). Spin up is shown along the positive, spin down along 
the negative axis. The black line is the QDOS of the right ferromagnet, treated within the NSDA, 
and the red broken line shows the QDOS of the left ferromagnet in mean field. 

shift is mainly influenced by three quantities: the strength of the hybridization itself {e^p 
and eMi in our case), the energetic distance between the two bands (the closer they are to 
each other, the stronger they will be repelled) and their spectral weight (higher spectral 
weight leads to stronger repulsion). In the upper part of Fig. [3] we plotted a typical QDOS 
for the NSDA without applied voltage for parallel alignment of the two magnetizations. 
The dashed line represents the density of states of the left ferromagnet. The splitting of the 
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right QDOS into lower and upper Hubbard band at ~ and E ^ Ur is clearly visible. 
Additionally there are contributions of the insulator at E Tq j = 5 eV which are due to 
the hybridization. What happens when a voltage is turned on depends on its sign. The left 
ferromagnet FMl is not influenced by the voltage, so its QDOS will be the same. Let us 
first discuss the case V > where both spin bands of the right ferromagnet FM2 are shifted 
to lower energies. But due to the repulsion between the spin up bands of FMl and FM2 the 
shift of the spin up band will actually be stronger than for spin down, where the repulsion 
is much weaker. So a positive voltage leads to a stabilization of the magnetization. This 
effect is enhanced by the current. For positive voltages it will flow from left to right. Since 
the left ferromagnet is fully polarized there are only spin up electrons tunneling into the 
right ferromagnet. These are the reasons for the slight increase of magnetization in Fig. [2] 
between points A and H and also an explanation why there can be no switching in this case. 
On the other hand, if we apply a negative voltage, < 0, the right spin bands will be shifted 
to higher energies. For the same reasons as discussed above the shift of the spin up band 
will be enhanced by the hydridization. Thus the difference between the centers of gravity 
of both spin bands will be decreased. If the hybridization strength is sufficiently large this 
additional shift together with the selfconsistency will be enough to push the spin up band 
above the spin down band and thus the magnetization changes sign. The selfconsistency is 
important since it will enhance the shift, because the occupation in one band depends on the 
occupation in the other one. In this case current flows from the right to the left lead. Since 
there are only spin up states available for tunneling in the left ferromagnet, the tunneling 
current flowing out of the switching magnet will consist of spin up electrons only. This leads 
to an additional decrease of the right magnetization. This explains the behaviour shown by 
the parallel aligned curve in Fig. [2l 

The antiparallel case can be explained in a very similar way. In the lower part of Fig. [3] the 
corresponding quasiparticle density of states is shown, again without applied voltage. The 
obvious difference to the case discussed above is that the center of gravity of the lower spin 
down band is below the lower spin up band. This is the reason for the negative magnetization 
of course, but it is also responsible for the reversed behaviour with respect to the applied 
voltage. In this case a negative voltage cannot push the center of gravity of the spin up band 
below the spin down band. It rather has the opposite effect because the repulsion pushes the 
spin up band to even higher energies and thus leads to a more stable magnetization which 
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can also be seen in Fig. [2] between points E and D. For a positive voltage the spin up band 
of FM2 moves below the spin up band of FMl so that the hybridization will shift it to lower 
energies compared to the spin down band. Again, if the hybridization strength is larger than 
a critical value this additional shift will be enough to reverse the two spin directions such 
that the magnetization changes sign. For the same reasons as discussed for the parallel case, 
a positive voltage will increase the magnetization, while a negative voltage has the opposite 
effect. Therefore the behaviour of the magnetization in Fig. [2] can be understood in terms 
of the quasiparticle density of states. 

In order to prove the explanation based on the hybridization we plotted a second magneti- 
zation curve in Fig. [2] with smaller hybridization strength eMi = 0.2 eV. Obviously in this 
case no switching occurs. There is only a slight change of magnetization. Starting from par- 
allel (antiparallel) alignment the magnetization is reduced for negative (positive) voltages. 
This is in agreement with the explanation given above. Since the hybridization is weaker 
the repulsion between the bands is also reduced. It is not strong enough to push the spin 
down band above the spin up band or vice versa. Thus the direction of magnetization is not 
changed. The current density through the junction is closely linked to the coupling strength 
between the materials^: smaller eMi corresponds to a weaker current. From the results 
shown in Fig. [2] we can conclude that in order to switch the magnetization the current has 
to exceed a certain value. 

To summarize, we presented a selfconsistent calculation of the voltage-dependent magneti- 
zation in a magnetic tunnel junction within a microscopic non-equilibrium framework. The 
magnetization shows a hysteresis behaviour similar to that seen in experiments. The reason 
for this effect was explained to be the hybridization between left and right ferromagnets 
which could be seen with the help of the quasiparticle density of states. It should be noted 
that the behaviour discussed above does only appear for very special parameter sets (such as 
low band occupation, small U) when one uses the mean field approximation for the right fer- 
romagnet. This seems reasonable because it is known that mean field strongly overestimates 
the stability of ferromagnetism. Thus it should be more difficult to switch the direction of 
magnetization. We have to conclude that higher correlations seem to be an important factor 
when describing current-induced switching of magnetization within this model. One might 
argue that the Kondo peak is missing in the NSDA which should have considerable influence 
on the magnetization. However, we investigated the strong-coupling regime only, where it 
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is known that the Kondo peak does not play a major role. On the other hand it would be 
a very interesting expansion of the model to examine its weak-coupling behaviour. Another 
important extension would be the inclusion of spin-orbit coupling which is widely believed 
to be the microscopic origin of phenomenological damping effects^- which play a crucial role 
in the macroscopic description of switching of magnetization. 



APPENDIX A: NON-EQUILIBRIUM SPECTRAL DENSITY APPROACH 

The basic idea behind the nonequilibrium spectral density approach (NSDA) is to choose 
the selfenergy in such a way, that the first four spectral moments 

oo 

Mit = -l I E-l^Gl^^{E) (Al) 

— oo 

of the spectral density are reproduced exactly. The moments are calculated with the help 
of the following exact relation: 



I = ^ E ([[... [c.,., H]...., HU [H... [H, c+ J_ . . . (A2) 



where the total number of commutators on the right hand side must be equal to n. Inserting 
the Hamiltonian ([1]) into this expression yields after some calculation: 

Mlt = 1 (A3) 

= + UrU-^ (A4) 

= 4, + 2f/Rek,n_. + f/>_. + e^j + tip (A5) 

(3) 



+[/||(2 + n_^)ek^n_, + - + Uln^, (A6) 

The moments of the transport selfenergy AJ^^^(i?) can be derived in the same way. One 
gets: 

D^l = (AT) 
= 4// + 4p (A8) 
D^l = el.Toj + elpTo,p (A9) 
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The band correction is given by 

+ ^ E E(^^x^«<<x-<.c.«-.(2n.,.-l)) (AlO) 

X=I,PiRix 

where Tj^j^ is the hopping integral between lattice sites Rj^ and Rj^. The two higher 
correlation functions can be reduced to single-particle lesser Green's functions^^. We find 

CO 

■GU.-AE) + A<^^AE)GU.-^(E)] (All) 

and 

oo 

{ct-.c,,^.<.) = -^^E^^'"^''^""'''"^ / dEJ:i,^^AE)G<^^AE). (A12) 

-oo 

The non-diagonal lesser Green's function G^^^^^{E) = i{d^^^Ck,^„) is closely related to the 
right Green's function and the transport selfenergy: 

E E GU.^exn = GL,.A<^, + G<^.A^^, (A13) 

X=I,P kx 

Putting all these expressions into the band correction leads to the result in Eq. [HI 
The Dyson equation of the right ferromagnet reads 

EGl^^iE) = 1 + (6k, + J^l^AE) + K.AE))Gl^^{E) (A14) 

Inserting the high-energy expansion for both selfenergies and the Green's function 

oo 

GIAE) = Y,^ (A15) 

n=0 
oo jj{m) 



oo Q{rn) 

^lu^^E) = E 

m=0 



yields a system of equations for the unknown moments G^^ of the interaction selfenergy. 
It can be solved by sorting according to the order of 1/i? and the use of the moments given 
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earlier in this appendix. The results are quite simple: 



= Uru^'^ (A18) 
C^l = f/^n_.(l-n_.) (A19) 
Cki = U^ri.^il - n_,)B_, + f/|n„.(l - n^^f (A20) 



These expressions are formally identical to the equilibrium case, therefore the selfenergy will 
also have the same form Eq. ([7])^^. For high energies it is acceptable to neglect higher order 
terms of the expansion in Eq. (lA17p . Thus: 

(1) 



^(0) C'k^a 



kjjcr ^ „(2) 



TT E — Tqh — B_„ /^o1^ 
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